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ABSTRACT 

The time dependent tidal interaction of the Local Group Dwarf Spheroidal 
(dS) Galaxies with the Milky Way (MW) can fundamentally affect their dynam- 
ical properties. The model developed here extends earlier numerical descriptions 
of dS-MW tidal interactions. We explore the dynamical evolution of dS systems 
in circular or elliptical MW orbits in the framework of a parametric oscillator. An 
analytic model is developed and compared with more general numerical solutions 
and N-body simulation experiments. 

Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics 

1. Introduction 

The Local Group dwarf spheroidal (dS) galaxies may be sensitive probes of the chemical 
and dynamical formation history of the Milky Way (cf Harbeck et al. 2001). They have 
also been used to constrain dark matter models and early structure formation scenarios in 
the Universe (cf Lokas 2002). 
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Controversy over whether the dS are equilibrium dynamical systems (cf Mateo 1998) or 
tidally driven non-equilibrium objects (Kuhn and Miller 1989, Kroupa 1997) persists. The 
primary observational evidence for tidally dominated dynamics comes not only from the 
dS high velocity dispersions and highly elliptical shapes but also from measurements of a 
spatially variable velocity dispersion (Kleyna et al. 2002) and significant stellar populations 
beyond the dS nominal tidal radii (Kuhn et al. 1996, Martmez-Delgado et al. 2001). While 
there is undisputed evidence that the dS in Sagitarius is tidally disrupted (Ibata et al. 1994) 
some practitioneers have devised more complex multiparametric models which could formally 
account for the kinematic data in some of the high M/L dwarfs using dark matter (Kleyna 
et al. 2002, Wilkinson et al. 2002, Lokas 2002). 

A theme of the tidal excitation argument presented here is that the gravitational interac- 
tion scale for transfering energy to the dS stellar system from the dS-MW orbit can be much 
larger than the dS classical tidal radius. In those systems where the stellar crossing time can 
be comparable to the "period" of the external tide, the simple static arguments for equilib- 
rium break down. Early calculations (Kuhn and Miller 1989 - henceforth KM) demonstrated 
this possibility using particle-mesh gravitational N-body simulations in idealized circular dS 
orbits. 

This paper develops an analytic formulation of the dS-MW tidal interaction problem. 
Once again, by treating the response of the dS in terms of galaxy oscillations, but now 
incorporating the tidal interaction through the Mathieu equation, it is possible to understand 
the evolution of these tidally interacting systems in non-circular orbits far from the simple 
resonance condition. As a parametric oscillator we find that a dS stellar system can be 
"inflated" even when its characteristic pulsational spectrum is not tuned to its MW orbital 
dynamics. 



2. Parametric tidal interactions 

An oscillator driven by a periodic forcing function can sometimes be described by the 
Mathieu equation 

x + uj 2 (1 - £ cos2cc;t) x = (1) 

Here x(t) describes the oscillator amplitude function, ojq is the oscillator resonant frequency, 
uj is the driving frequency and e is a small constant. Solutions to equation (1) are most 
easily developed by a perturbative multi-timescale analysis (Bender and Orzag 1978). One 
finds that, in general, growing (unstable) oscillatory solutions exist for uj = u /n, n G N. 
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An important property of equation (1) is that even for a weakly damped (high 'Q') 
oscillator a significant, non-zero, frequency range for the driving force leads to growing 
unstable solutions. For small intrinsic damping this result is independent of the frequency 
width of the natural resonance. For example, near the first resonance u between ou (1 — e/4) 
and a>o (1 + e/4) leads to instability. The second resonance near uj = ujq/2 is unstable for a 
narrower frequency range 

— e 2 u > u-u /2 > ^e 2 u (2) 

We show here that the internal dynamics of stellar systems in the MW's gravitational 
environment can be described as a two-dimensional coupled parametric oscillator. These 
equations can be solved analytically using multi-scale perturbation methods to explain the 
dynamical and morphological properties of the Local Group dS. 



2.1. Tide geometry 

Each of the dS orbit the MW at a distance which is long compared to the galactic disk 
scale length so that it is reasonable to treat the MW as a purely spherical central force 
problem. We define the MW central force 

F MW =F(V) = F(r)r (3) 

Letting ~r describe the displacement between the MW and dS centers and 8r be the dis- 
placement to a dS star from the dS center we write the MW tidal force on the star at lowest 
order as 



AF MW = F (r+5rj - F ( r ) 

= F (||~^ + <Sr||) - F(||7*||)] r 
+ sin (Ay?) F( || 71) £ 



(4) 



Here A(p is the central angle between a dS star, the MW, and the dS center and f and tp are 
unit vectors in the radial and central angle directions (Figure 1). Since the spatial extent 
of the dS is small compared to its distance from the galactic center we expand equation (4) 
using 
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(||^ + ^||) -^(ll^ll) 
sin (Aip) 



— — (r) or ■ r 
dr 

- 5r • (p 
r 



(5a) 
(5b) 



We define a cartesian (xy) coordinate system centered on the dS which contains the MW 
and the star so that we can express 5r, r and (p in term of x and y. 



5r = xx + yy 
r = cos <p x + sin tp y 
(p = — sin ip x + cos ip y 



(6a) 
(6b) 
(6c) 



The MW orbital time dependence is contained in r and tp. To this order of approximation 
these equations also describe the three-dimensional stellar problem where (x, y) are the 
projected coordinates in the dS-MW orbital plane. The tidal force perpendicular to the 
orbital plane yields harmonic motion which is not coupled to the MW orbital motion and is 
ignorable. Thus the tidal force becomes, 



dF 

dr 
dF 

dr 



- F ( r ) ( 2 
r) [ x cos^ ip + ^ sin 2ip ) H — ism ip 



y 

2 



(r) ( : 



y sin 2 tp + — sin 



2y?) 



r 

F(r) 



H — - (i) cos 2 tp — ^ sin 2ip j 



sin 2<p 



(7a) 
(7b) 



The form of the MW potential determines ^- and F. For a Keplerian potential, it is 
dF _ _2 E(r) anc j f or logarithmic case it is 4r- = — Hence we can write a simplified 

dr r ° dr r ^ 

form of the tide for these two potential forms: 
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dF 

dr 



cos 2ip — — + y sin 2ip j x 
+ ^— y cos 2(/9 — ^ + x sin 2</? j y 
a; cos 2ip + y sin 2<^) a; 

+ (— y cos 2(/9 + x sin 2tp) y 



(8a) 



(8b) 



The numerical factors are of course different for different potential assumptions but the form 
of the resonant solution remains unchanged. It seems likely that the MW potential is actually 
logarithmic in the domain of most dS so that we assume the form 



= v c In — 

Jo 



(9) 



2.2. Solution of the coupled Mathieu equations at first order 

2.2.1. Equation of motion for a star near the core 

In the dS reference frame a star experiences forces due to the MW tide and the dS 
potential. The dS potential generates a radial acceleration, a r = —GM(r)/r 2 , where M(r) 
is the mass enclosed within a spherical volume of radius r centered on the dS. Near the 
center of the dS we approximate M(r) p 4:Trr 3 /3 so that we can write a r ~ — u 2 r with 
Ld 2 = AnGp /3. In this case the total dS-frame acceleration of a star can be written as 



dt 2 m 

With additional definitions, 



d 2 7> 2 _> AF 

-cu 2 r + (10) 



k = 



e = 



1 dF 

m dr 
k 



(11a) 
(lib) 
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where the circular frequency of the dS around the Milky Way is uj c = v c /r. Combining these 
terms we finally obtain a recognizable set of coupled equations that describe the behavior of 
a star in the center of a dS as projected in the (xy) plane. 



Notice that the coupling between x and y motion results from the anisotropic form of 
the galactic tide where the expansive and compressive tidal contributions are both important 
in determining the dS dynamical response. 



Individual stars also experience forces due to the dS gravitational perturbations but 
equations (12) neglect temporal variations (due to p(t)) in the dS potential. Thus there is 
an additional acceleration term which should be included here at late times in the dS-tide 
evolution when the temporal variability of the dS density is important. 

The resulting oscillations in the mean stellar orbits are not included above but can be 
described in a manner which is analogous to polytropic stellar normal mode oscillations. 
This description of galaxy dynamics in terms of linear oscillations was the basis of KMs 
prediction that Local Group dS could be tidally inflated even when the stars were contained 
within the static tidal radius. Miller and Smith (1994) and VanderVoort (1999) have now 
demonstrated that weakly damped or even growing finite amplitude normal mode galaxy 
oscillations are a physically interesting feature of realistic galaxy dynamics simulations. 

The lowest order oscillation mode should have no radial nodes - as is the case for radial 
oscillations in stellar polytropes (cf Cox 1980). VanderVoort (1999) showed that the total 
Lagrangian radial displacement summed over all stars in an N-body system exhibits a very 
similar oscillatory behavior. This is not entirely intuitive given that the collisionless orbital 
motion of stars in a galaxy occurs over the length scale of the galaxy, while a mass element 
in a stellar polytrope oscillates only over the small scale of the oscillation displacement 
amplitude. Nevertheless VanderVoort proves that the fundamental N-body radial eigenmode 
for the mean Lagrangian radial particle displacement, £ , satisfies 




(12a) 
(12b) 



2.2.2. Galaxy oscillations 
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dt 2 



= -<V £ (13) 



with lu 2 = —W/I = G J M (r) p(r) r dr / Jp(r)r 4 dr and £ oc ~r. Here W and / are the 
potential energy and moment of inertia, while p(r) and M(r) are the galaxy density and 
enclosed mass distribution and G is the gravitational constant. Realistic density distributions 
lead to frequencies that are of order n^Gpo where po is the central density. 



2.2.3. Circular orbit 

If we launch the dS on a circular orbit around the MW, then ip = ut and equations (12) 
describe a set of two-dimensional Mathieu equations. We now show that these equations 
exhibit resonant instability around u = ujq. To obtain a valid solution near the resonance, 
we use a multi-timescale analysis in terms of t and t — et and develop the relevant variables 
in terms of the perturbation parameter e. 



u>o = u + eui (14a) 
x(t,r) = x (t,T)+ex 1 (t) (14b) 
y(t,r) = y {t,T)+e yi (t) (14c) 



The equations obtained from the lowest order (e°) for x and y Q are simple harmonic 
oscillators with frequency u as we would expect in a harmonic system. From the lowest 
order xq and yo solutions and first order equation we can construct a refined solution. Using 
complex notation we write 



x (t,r) = A(t) e i0Jt + c.c. (15a) 
y (t,r) = B(r)e^ + c.c. (15b) 



Gathering the terms of order e 1 , we obtain differential equations for x\ and y\ 
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d 2 xi o 
d 2 y 1 



—2ujijj\ x + cu 2 cos (2 cut) x + u 2 sin (2 cut) y — 2 



dt 2 



+ cu 2 y 1 = —1ujijj\ y — J 2 cos (2 cut) y + cj 2 sin (2 cjt) x — 2 



d 2 x 
dtdr 
d 2 y 
dtdr 



(16a) 
(16b) 



These are combined with the lowest order solutions for xq and y . We construct a solution 
for A(r) and -B(t) using the well-known ansatz (Fredholm alternative) of forcing all har- 
monic terms which lead to unstable (secular) solutions to vanish. Obtaining a solution with 
unbounded coefficients implies that there is no stable oscillatory solution for x and y. Thus 
we force the coefficients of e lut and e~ lu)t on the RHS of equations (16a) and (16b) to vanish 
to obtain two differential equations (and their conjugates) for the functions A(t) and B{r): 



-2. 1 X(r) + |x- W + ^B-(r)-2i^ = 



(17a) 
(17b) 



We solve these four coupled differential equations in terms of the real and imaginary 
parts of A and B, i.e. (ft (A) , 3 (A) , ft (B) , S (B)). This yields the matrix equation 



where, 







dX 
d7 


= MX 
(ft (A)\ 












X = 


91(A) 
ft(B) 










( ° 


- (4wi + cu) 


—cu 
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(4cui — cu) 










cu 




4 


—cu 










-(4wi 






\ 


cu 


(4wi + cu) 





/ 



(18) 



(19) 



(20) 



This equation is solved most readily by finding the eigenvectors of M. Taking 



H 



cu + 2cu\ 
cu — 2cu\ 



(21) 
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the eigenvectors of M of respective eigenvalues 



i/ai 2 -4wi 2 \/w 2 —Au) 1 2 . i 

g ) ~ 2 ' 1CJ 1 aI1C ^ _ 1CJ 1 are 



1 
1 



( H\ 

1 
1 



— i 
i 

V 1 / 



, and 



i 

-i 

V 1 / 



(22) 



The solution to equation (18) is obtained in terms of four complex constants of integration 
C = (a, (3, 7, 5) and the eigenvectors as: 



a exp 
(5 exp 



\Jtjj 1 - Au{< 



y/u 2 - AuJi 2 



-H^(A) + Q(A) + ^(B) + HQ(B) 



7 e^ 1 ' 



(23a) 



= H$t(A) + S(A) + -H S(B) (23b) 

U(A)-i^(A)+m(B) + Q(B) (23c) 
M(A)+i%(A)-m(B) + Cs(B) (23d) 



At late times the growing exponential term dominates, and these equations simplify to 
yield expressions for A and B, 



A{r) = 
B(t) = 



a 

Ih 

a 

Ih 



exp 
exp 



y/u 2 - 4cJi 2 



yJuS 2 - AuJi 2 



et 



et 



(-1 + iH) 
(H + i) 



(24a) 
(24b) 



Evidently A and B define the mean position of a star in the dS 

' \Juj 2 - Auoi 2 



x(t) = a' exp 
y(t) = a' exp 



\/uo 2 - 4cji s 



et 



et 



cos (ut + 0) 
sin (ut + 0) 



(25a) 
(25b) 



Equations (25a) and (25b) show that instability occurs if u > \2ui\. Thus if the circular 
frequency of the orbit satisfies 
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ujq — u c \ < u) c e/2 



(26) 



the dS will be tidally excited. The constants C determine the initial position and velocity 
of a given star in the dS. Notice that the growing eigenmode solution involves stars oriented 
along a line that makes an angle 9 with respect to the x axis where tan# = —H. Thus, we 
find in general that stars expand away from the center of the dS depending on the difference 
between the dS orbital frequency and c^o, and on the perpendicular distance a star makes to 
a line in the (xy) plane of the dS at angle 9 with respect to the line-of-sight toward the MW. 

Extending this calculation to second order in e reveals that, unlike the simple one- 
dimensional case, the next resonance near uq/2 is stable with no secular amplitude variation. 



Equations (25a) and (25b) assume a circular orbit. Their applicability is also restricted, 
depending upon the central concentration and oscillation frequency u of the dS. For exam- 
ple, the perturbation parameter, e = uj c 2 /ujq 2 , can approach unity for a circular orbit in a 
logarithmic potential. Under these conditions the static tide force alone is strong enough to 
disrupt the dS. Since the dS are generally not in circular orbits it is important to extend 
the model to elliptical orbits. In a logarithmic MW potential we can recover the form of the 
Mathieu equation with the epicycle approximation. We define an "eccentricity" e from r(t) 
where k = \/2uj is the epicycle frequency for a logarithmic potential. It follows that (Binney 
and Tremaine 1987), 



2.24. Elliptical orbit 



r(t) 
<p(t) 



r c (1 — e cos (nt)) 
uj c t + e \/2 sin (nt) 



(27a) 
(27b) 



Here k = (v c /r(t)) 2 pa uj 2 (1 + 2e cos (nt)) where v c = r c co c is the circular velocity of the 
potential. We note that with this convention for r(t), the dS is at perigalacticon at t — 0. 
Replacing k(t) = e(t)uo 2 and <p(t) in the equations (12) gives, at first order in e, 
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x + ujq 2 ^1 — e c cos (2cu c t) — b ± e e c cos (2a± u c t)^j x = 

coo 2 ^e c sin (2cu c t) + ^ 6± e e c sin (2a± c; c t)^ y (28a) 
i/ + c^o 2 ^1 + £ c cos (2u; c i) + b ± e e c cos (2a± oo c t)^j y = 

u 2 ^e c sin {2u c t) + 6± e e c sin (2a± w c t)^ x (28b) 

where a± = 1 ± \/2/2, 6± = 1 ± y 7 ^ and e c = u c 2 /uj 2 . Equations (28) have a form identical 
to equations (12) where uj is replaced by a ± uj c and e becomes b± ee c . Evidently two sets of 
resonant frequencies occur when uj c {1 + V2/2) or co c (l — V2/2) is near cu . These frequencies 
are sufficiently different that we treat them as independent resonant solutions. It follows 
that the form of the growing mode late-time solution for x(t) and y{t) in an elliptical orbit 
is 



x{t) = a' exp 



y(t) 



a exp 



a ± u) c ) — 4cjx 2 



a ± uj c ) 2 — 4cJ]/ 



\b±\ ee c t 



ee c t 



cos (a±u c t + 9±) 



sin (a±u c t + 9±) 



(29a) 



(29b) 



where lo\ has been defined by u = a±u c + b±ee c uj 1 . Here tan6* + = — H + and tan#_ = l/i7_ 
where H ± has the form of equation (21) with u replaced by a ± u c . Comparing to the circular 

2 

orbit solution we see that the expansion parameter e becomes e b ± so that instability 

2 

occurs when \u — a±u c \ <u ^e \b ± \/2. This expression is comparable to the circular 
orbit solution (26) but resonance is "easier" to reach even as the strength of the tidal force 
compared to the binding force (which is given by e c = k c /uj 2 = (uj c /ujo) 2 ~ l/a ± 2 ) can be 
weaker by a factor of 3 or more compared to the circular orbit. 

Expanding (p(t) and e(t) to order 3 in e yields additional resonant behavior. Analogous 
to the above expansion we find two new sets of excitation frequencies. The additional terms 
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in the differential equation for x(t) and y(t) above have the form d,2±e 2 e c cos (2b±u c t) and 
d 3 ±e 3 e c cos [2c ± uj c t) - with corresponding sin(-) terms as in (28). Here c± = 1 ± | \/2 and 
d n ±, n G {2,3}, are numerical factors of order unity. Thus, the higher order terms lead to 
excitations at frequencies where Uq is b± uj c (order 2) and at c ± u c (order 3). The growth 
rates of these modes and instability range for ujq are smaller by corresponding powers of e 
so we expect to excite higher order modes only as the orbital eccentricity increases. 



2.2.5. Projection on the line of sight 

An interesting, and perhaps observable signature from the simulations, is the elliptical 
shape caused by the strong geometric asymmetry in the growing mode solution for the 
distance of a star away from the center of the dS. Here we consider how the induced elliptical 
shape varies with respect to the line-of-sight from the MW center-of-force toward the dS 
center. Consider, for the sake of generality, that the excitation is encountered at a frequency 
au) c where a can be chosen from the set {1, a±, b±, c±}, with the corresponding exponential 
term E a (t). Thanks to equations (6b) and (6c), we obtain (where t = corresponds to 
perigalacticon) 



(xx + yy)-r = E a (t) cos ^(a — 1) uo c t + (30a) 
(xx + yy)-(p = E a (t) sin ((a - l)u c t + 6») (30b) 



In the circular case (a = 1) the dS is expanded along an axis rotated away from the 
line-of-sight to the MW by an angle 9 with tan# = —H. For elliptical orbits the bar appears 
to rotate. At order n the bar formed by the tidal interaction will be turning at a frequency 
nn/2. Thus as the dS orbits between successive perigalacticon the bar will turn | times. 



3. Numerical Calculations and Experiments 

These analytic calculations suggest that a broad range of dS-MW orbit dynamic con- 
ditions can lead to parametric tidal dS excitation. Evidently this excitation can occur even 
when the orbital and dS galaxy resonant frequencies are widely "mistuned." Although our 
analytic calculations simplified the force equation terms from the exact spatial and tempo- 
ral dependence of the dS self-gravity, we can show by numerical methods how many of the 
properties are retained in more realistic self-gravitating stellar systems. 
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In the discussion below we assume a system of units where the gravitational constant 
G — 1. We fix one numeric length unit to a physical scale of 1 kpc and we've taken the 
dS to have a total initial mass of 10 6 solar masses. With the gravitational constant set to 
unity, this implies a numerical time unit of 415 Myr. The logarithmic external Milky Way 
potential has been chosen to have the form 0(r) = t> c 2 log(r /ro) where we assume a circular 
velocity v c = 220 km.s -1 (or 100 numeric units). 

We first consider the problem of a more realistic non-harmonic, but still static, dS 
potential/density model. This is addressed by numerical integration of the modified two- 
dimensional Mathieu differential equation. Finally we consider direct N-body calculations 
to explore the collective and time-dependent effects of the dS potential as it responds to an 
external tide. 



3.1. Plummer Potential Solutions 

The Plummer model, 0(r) = —Mj \J r 2 + b 2 (cf Binney and Tremaine 1987), is a rough 
but useful approximation for a self-gravitating dS stellar distribution. We use this model 
both to initialize our numerical simulations and in the differential equation integrations 
below. For example, in this case the local cartesian x acceleration of a star can be written 
as x, = —xujq 2 ^ 2 ^ 2 ^ /2 where u 2 = 4np /3 and p is the central density. Notice that the 
parameter b conveniently describes potentials which range from completely unbound (6 = 0), 
to homogenous, purely harmonic {b — > oo) models. Substituting for the harmonic dS force 
term in eqs. (12) we obtain 



b 3 

9 9 9 / \ 

x + xujo 3/2 — xsujq cos2(/9 = yeuo sin2y9 (31a) 

{b 2 + r 2 )' 1 

b 3 



y + yuiQ 2 rrr + y e u 2 cos 2ip = xeu 2 sm2(p (31b) 

(b 2 + r 2 ) 1 



We were unable to construct analytic solutions to equations (31), except for the cases 
6 = 0, or b — > oo. Notice that this more general equation, preserves the form of the two- 
dimensional parametric oscillator but with a spatially variable "resonance" frequency w p 2 (r) = 
u 2 b 3 /(b 2 + r 2 ) 3 ^ 2 and with e — > eu 2 /lj p 2 . For example, stars on circular orbits with ini- 
tially constant r should be tidally excited at the first Mathieu equation resonance condition 
when the dS orbit satisfies uj = ui p (r). We can expect that, to the extent that individual 
stellar orbits aren't circular and there is a broad range in u> p (r), the dS tidal response should 
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be a broad function of its orbital frequency, uj. For finite b we also note that u> p (r) < ujq 
so that resonant behavior must occur at lower frequencies than the dS harmonic oscillation 
frequency. 

We can explore this solution space numerically. To distinguish stable and unstable 
solutions as a function of orbital frequency uj we numerically integrate eqs. (31) using a 
fourth-order Runge-Kutta algorithm. For example, according to the analytic solution, with 
b — > oo, Uq = 4, and e = 0.1 we should obtain solutions for r(t) = x 2 + y 2 which are 
growing only for dS orbits with circular frequencies uj ~ ujq. To verify this we integrate an 
ensemble of solutions with uj between and 8 over t = to 50 and compute the mean square 
radius over time for each solution. Here x(0) = 1 and the stellar orbits were chosen with a 
zero initial velocity condition to yield radial orbits. Figure 2 plots this mean solution radius 
for each dS orbital frequency. 

Figure 2 illustrates several features of the analytic parametric solution. First, instability 
occurs when the orbital frequency is close to the dS resonant frequency (u p ) and, unlike in the 
one-dimensional parametric oscillator, the next resonance where the orbital frequency is one- 
half of the oscillation frequency is stable. This agrees with our analytic result. It is perhaps 
surprising that a circular orbit with half the orbital frequency of the dS resonance does not 
lead to instability - since the tide force has 180-degree symmetry. This is a consequence of 
the compressive and expansive parts of the tensor tide interaction. 

The instability frequency width is determined by the tidal strength and form of the MW 
potential (through e in eq. (26)). With uj = uj c = 4 and e = 0.1 we expect a resonance 
width of 0.4 - which is confirmed in figure 2. The growth rate of r(t), as determined by eq. 
(25), is also reproduced in the numerical calculation. For uj± = we expect a stellar orbit to 
expand like expuj et/2. = exp0.2t which by t = 50 yields the amplitude plotted here. 

The other extreme, where 6 = 0, corresponds to a dS which is completely unbound, i.e. 
where the dS self-gravity is unimportant. For example, at late times a dS might be tidally 
distended so that its self-gravity becomes negligible. A naive interpretation of the MW effect 
here might be to conclude that unbound dS member stars should be tidally accelerated to 
large distances from the dS center-of-mass. Figure 3 shows the logarithm of the rms stellar 
orbit radius for the unbound (b = 0) case with the same initial conditions, range of dS-MW 
orbital frequencies, and integration parameters as in fig. 2. We see that for static and slowly 
rotating tides the stellar system does indeed rapidly blow up - saturating our numerical 
dynamic range before t — 50 for orbits near uj = 0. Analytic solutions using computer aided 
symbolic manipulation imply that the leading growing mode term in the solution varies as 
exp [a/ eoj 2 — uj 2 t\ . Interestingly, as this shows, for dS orbits with short enough periods 
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satisfying 

uj > \/suj = Vk (32) 
localized near the dS center-of-mass, despite being unbound to the dS. 

Intermediate potential models (0 < b < oo) yield unstable orbits depending on the 
initial conditions of the stellar orbit. The dS orbital frequency which leads to instability of a 
dS star depends on the initial radius of the star with respect to the dS center-of-mass in units 
of b, and the dS central density which determines ujq. The previous two cases have shown 
that dS stars which begin their orbits at normalized distances of and oo develop secular 
instability for dS-M W orbits that satisfy, respectively, uj ~ uj and uj rs 0. Intermediate cases 
lead to intermediate resonant frequencies. 

The top panel of figure 4 plots the orbital expansion of a star starting from r/b = 1 using 
the same integration parameters (ujq = 4) as the previous Runge-Kutta solutions. Secular 
growth now occurs at lower frequencies near uj m 2.8. In a Plummer potential-density model 
the half-mass radius occurs at 1.46. The lower panel of figure 4 demonstrates for this radius 
how the dS resonant orbital frequency continues to decrease for stars at larger normalized 
dS distances. Evidently stars near the half-mass radius are tidally excited when uj — 1.1. 
The initial position of a star within the dS potential helps to determine its fate in the 
time-dependent MW tide. 

These solutions describe a system where the tide amplitude e and the orbital frequency 
uj are independent. A logarithmic MW potential satisfies euj 2 = uj 2 . Thus for increasing 
orbital frequency, e in eqs. (31) is also increasing. Qualitatively we find that this increases the 
orbital frequency domain over which resonance expands the dS. Systems which are resonant 
with uj rs uj p for fixed e in the above discussion are now tidally disrupted for uj > uj p . 



3.1.1. Elliptical Orbits 

Elliptical orbits are also described by eqs. (31) by taking e — > e(t) and using the actual 
form of if{t) and e(t) for an elliptical orbit. Here we assume a logarithmic form for the MW 
potential so that e and uj are not independent. 



tp(t) = l + ev^sin/tt (33a) 

= ^4 T, v7 (33b) 

v ; uj 2 (1 - ecosKt) 2 v ; 
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where k = a/2 uj is the epicyclic frequency. 

We compute the mean distance of a star from the dS center as we did in the previous 
section. We assume a constant perigalacticon distance with e = 1 — uj /2, take uj = 6.7, 
x(0)/b = 1.4, and let vary from 1.0 to 2.0 so that we can compare with the results of 
section 3.2.3. Figure 5 shows that a broad range of circular frequencies lead to parametric 
excitation and, as we will see, this general behavior also appears in more realistic self- 
gravitating systems. We also note that a circular orbit (uj = 2) can be less dispersive than 
an elliptical orbit (uj < 2) even if the tide is always stronger in the circular case. 

3.2. Direct N-body Calculations 

In order to account for the self-gravity and collective dynamics of many stars within 
the dS we have used a direct N-body calculation. The numerical simulations were computed 
using the modified treecode vl.4 (Barnes 1990). Typically between 1000 and 10,000 particles 
were introduced into simulations which included several different "external" forces designed 
to study variously dS internal oscillations and the tidal coupling of the stellar system to an 
imposed large-scale logarithmic external MW potential. 

The effects of numerical viscosity and small spatial scale potential fluctuations can be 
important in these calculations and we used a range of integration time step and potential 
"softening" parameter to understand this dependence. Most of the simulations were done 
with a time step of 0.008 units. We explored the effects of the interparticle potential softening 
parameter over values ranging from 0.025 to 0.3. Significant damping of oscillations appears 
to be minimal in these models with softening parameter a > 0.1. 

With a crossing time of about 0.2 time units the relaxation timescale for N = 1000 
simulations was at least 2 time units and much larger for softer potentials. Many of our 
results only depend on measuring differential changes between simulations, e.g. to determine 
the frequency dependence in the particle ejection rate between models with different orbital 
frequencies. Some results are computed from relatively long model simulation runs. 

Several different initial particle configurations were generated although most of the 
results we present here used the treecode Plummer model realization (cf Aarseth et al. 1974 
as implemented by Barnes 1990). The rms particle radius was typically 0.1 kpc and the 
Plummer parameter was typically b = 0.16 for a softening parameter of a = 0.1. To achieve 
this an initially isotropic velocity dispersion of about 1 km.s" 1 was used and allowed to relax 
before imposing any external time-dependent tidal forces. 
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We have not tried to reproduce any of the dS in detail since their MW orbital uncer- 
tainties and the non-uniqueness of the model solutions make detailed comparisons difficult 
to interpret. Nevertheless we believe these N-body calculations are most representative of a 
dS like Draco or Ursa Minor. 



The parametric galaxy model described by eqs. (12) neglects the effect of the oscillating 
potential of the dS. Thus, we can expect normal mode galaxy oscillations to dominate even 
the growing parametric tidal modes if the system's normal mode frequencies are ever excited 
by the MW orbit. To explore this numerically we need to understand the normal mode 
oscillation spectrum of our dS system. 

Based on previous experience with particle-mesh calculations (KM) involving 10 5 parti- 
cles or more we expected free galaxy oscillations to be readily generated by small deviations 
from equilibrium in the initial conditions of the particle configuration. For example, Miller 
and Smith (1999) had difficulty damping and suppressing such oscillations. Even though the 
total energy in our simulations of isolated dS was conserved to within 0.03 %, we initially had 
difficulty detecting oscillations. To excite an oscillation spectrum we applied an impulsive 
radial stretching acceleration at t — and looked for the resulting ringing. From a tempo- 
ral fourier analysis of projections of the particle motion onto the fundamental eigenmode, 
£ (r) = ~r, we measure the oscillation spectrum. 

In detail we compute a timeseries by summing over particles (labeled with i) 



The logarithm of the temporal fourier transforms of Cj from two different simulation runs 
are indicated in Fig. 6. For an homogeneous sphere, the power spectrum should peak at 
the oscillation frequency ujq. Increasing the softening parameter both changes the peak 
frequency in the power spectrum and the peak amplitude response to the excitation. The 
virial equation (13) yields a reasonable estimate of the observed oscillation peak frequency. 
For example with softening parameter of a = 0.1 and using the observed N-body density 
distribution we compute a frequency of u = 3.54, which is consistent with the peak in Fig. 6 
(lower panel). 

The shorter timestep minimizes the integration errors and the intrinsic numerical damp- 
ing noise (Figure 7). Similarly by increasing the softening parameter we can decrease the 



3.2.1. Oscillations 




(34) 
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effects of integration error from the short range particle interactions. Figure 6 confirms 
that as the softening parameter increases there is an enhancement in the fundamental mode 
amplitude. 

It is difficult to compute the true physical damping of galaxy oscillations but comparing 
our results with Miller and Smith (1999) suggests that the direct N-body calculations have 
a larger numerical damping than particle-cell approximations. We expect the dS instability 
growth rates from these N-body simulations to be underestimated compared to the analytic 
and Runge-Kutta solutions. 

3.2.2. Parametric resonance 

In order to demonstrate Eqs. (25) for a dS in a circular orbit, we decouple the strength 
of the tidal force from its time dependence. We isolate the dependence of the tide amplitude 
and frequency by fixing the amplitude while independently setting the orbital frequency {uj). 
This is analogous to the approach KM used in their particle-mesh calculations. It allows the 
frequency dependence of the resonance to be distinguished from the effect of a growing static 
tide strength as the dS galactocentric distance is decreased in order to increase uj. Evidence 
of this resonant behavior is illustrated in Figure 8 which shows a snapshot of the particle 
positions for three different models corresponding to uj = uj p /2, uj p , and 2uj p . The direction 
of the instantaneous tide is indicated by the line on each plot. 

Figure 9 shows the effect of changing uj while keeping the tidal amplitude (determined 
by e) fixed. Here the total number of particles ejected from the dS is plotted versus time 
and frequency uj. Recall that uj p is approximately 2 numerical units and the tide amplitude 
corresponds to e = 0.5. Successive curves diplaced upward in this plot show how the ensemble 
of models have ejected more particles at later times. 

Several points are illustrated by fig. 9. The most important is that, as expected from the 
Mathieu equation solutions, the resonant behavior (measured by particle ejection from the 
dS) extends to frequencies relatively far from the effective resonant frequency uj p . According 
to our solution to Eq. (25) and for our choice of e we expect an instability frequency range 
Auj/uj p = e/2 = 1/4 which is consistent with the numerical results. It is also notable that 
there is no second order resonance at uj p /2, consistent with the analytic predictions of the 
coupled Mathieu equation model. 

We also find in the N-body solutions that for dS circular frequencies uj < 1.5 the ejected 
stars migrate toward r — > oo as the dS evolves. For uj > 1.6 the ejected stars (more distant 
than 1 kpc from the core) appear to remain in a bounded region of space, at least over the 
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duration of these simulations. This result is anticipated by eq. 32. In these simulations 
k = 2 so the cutoff frequency should be about u — 1.4, in fair agreement with the N-body 
calculation. 



3.2.3. Elliptical orbits 

More realistic orbital calculations also demonstrate how the dynamical tide interaction 
affects the internal dynamics of a dwarf spheroidal. We adapted the numerical calculations 
described above to describe a self-gravitating ensemble of masses orbiting within a logarith- 
mic external potential. We use a 1024 point dS galaxy characterized by oscillation frequency, 
uj 7 numeric units. Setting this dS onto a circular orbit at 50 kpc from the galactic center 
yields a stable system with 75% of the initial stars and with a particle loss rate of less than 
50 stars over 250 time units. 

To generate a family of elliptical orbits we launched the dS at 50 kpc along the x axis 
with variable speed vq > v c along the y axis in a constant galactic potential <f>(r) = v c 2 log(r). 
Thus in our units when uu c = 2 the orbit is circular and in general the eccentricity, e, varies 
like e = 1 — iv c /2. Thus, the perigalacticon of all runs is exactly 50 kpc and the mean tidal 
force on the dS is largest for the circular case (e = 0, uo c = 2) and decreases with uj c or as 
the eccentricity e increases. 

Figure 10 shows how the number of particles ejected from the dS varies with orbital 
ellipticity and time. Here the orbit time is measured from the injection of the equilibrated 
dS into the logarithmic potential. The broad frequency response of the dS in an elliptical 
orbit is confirmed from the Runge-Kutta solutions (Figure 5). It is interesting that elliptical 
orbits, with weaker tides than a circular orbit can lose more stars. Eccentric variable tides 
can eject stars where a stronger constant tide cannot. 

3.3. Velocity Dispersion 

The velocity dispersion of our parametrically excited dS galaxies varies by a large factor 
over their orbits and can be a strong function of the extent of the dS over which the dispersion 
is calculated. For example a general feature of these models is that the dispersion increases 
outward from the core of the dS. Nevertheless it is possible to inflate the dispersion by 
an order of magnitude or more depending on the dS resonant and orbital conditions. For 
example, in Fig. 11 we plot the instantaneous line-of-sight velocity dispersion for the model 
above with orbital eccentricity of 0.5. While this is a relatively high eccentricity orbit, it is 
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not particularly close to resonance (with a slow growth rate) and we plot the dispersion after 
the growing mode is established. It is clear that at late times the dispersion increases by more 
than a factor of 10 yielding virial mass estimates exaggerated by a factor of 100. Figure 12 
shows this model when the dispersion was 10 times larger than its initial configuration. This 
system is not particularly close to dissolution and it has lost only 300 stars out of 734 when 
it was launched in its orbit. It also seems to retain its central concentration and surface 
brightness profile. 

3.4. Interpreting galaxy morphology from N-body calculations 

Our analytic parametric oscillator calculations and numerical integrations of more re- 
alistic Plummer model systems showed that resonant excitation of dS stars should produce 
elliptical systems with their long axis rotated by an angle 6 away from the direction toward 
the force center in the plane of the orbit. For a circular or elliptical orbit this angle is cal- 
culated from tan# = — H (section 2.2.3). At the first parametric resonance this angle is 45 
degrees. Further from the resonance condition this angle varies from an orientation along 
the center-of-force direction to being perpendicular to it. The existence of a nonzero dS bar 
angle is evidence of parametric resonance. 

We found in section 2.2.5 that this angle also describes the dS elongation in an elliptical 
orbit at perigalacticon. During the orbit the bar of a dS should rotate with respect to the 
center-of-force direction with increasing frequency, depending on the order of the resonance 
in orbital ellipticity, e. 

We have computed the direction of the dS bar from the moment of inertia tensor of the 
dS stars. Figure 13 shows a series of snapshots from the simulation with e = 0.22. In this 
figure the center-of-force is always to the left and 8(t) is the angle between the two solid lines 
which intersect at the center of the dS. Figure 14 shows how the angle varies with time for 
this calculations. Each spike in the graph corresponds to one rotation of the bar. We find 
that the rotation rate increases with increasing eccentricity (as it should) as higher order 
modes in powers of e are excited. 

The bar rotation indicated in Fig. 14 is not strictly sinusoidal with the epicycle period. 
We also note that at perigalacticon our dS systems tend to have the bar directed toward the 
force center (6 = 0). Recall from our discussion of the Plummer model that not all stars in 
the dS are characterized by the same resonance frequency uj, so that stars near the center of 
the dS have larger frequencies than stars at larger distances. Thus for given orbital circular 
and epicyclic frequencies the response of many of the dS stars is non-resonant. Thus near 
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perigalacticon many of the stars (especially near the core of the dS) simply respond to the 
tide like any non-resonant fluid system. Our moment of inertia calculation is thus controlled 
by contributions from both resonant and non-resonant stellar components. This can lead to 
a more complex behavior of the bar pointing direction than our model predicts in detail. 

The outer parts of the dS (see Fig. 13) are curved as if H depends on distance from 
the dS. This is anticipated from the Plummer calculations since the effective u p decreases 
outward. The curvature in the bar at large distances from the dS raises the possibility that 
snapshot observations of real tidally excited stellar systems may provide useful constraints 
on the form of the dS potential. 

4. Comparing the parametric oscillator model with observations 

Our earlier comparisons of dS dynamics and morphology with a simpler resonant tidal 
excitation model depended on the commensurability of the dS oscillation frequency (and 
damping width) and its MW orbital driving frequency to generate elliptical systems with 
large velocity dispersions. The parametric oscillator model presented here more accurately 
describes the dS-MW interaction and is a more precise explanation for dS properties. In par- 
ticular their elliptical shape and sometimes large velocity dispersions is a natural consequence 
of this model. 

Parametric dS galaxy oscillations effectively increase the dS gravitational interaction 
cross-section with the MW. Properly accounting for the time- dependent tidal interaction of 
a dS with the MW, even when the dynamical frequency of the dS (cu ) is not commensu- 
rate with its orbital or epicyclic frequency around the MW, is critical to describing the dS 
dynamics. We have shown how the frequency domain near the dS fundamental frequency 
which describes resonance is not characterized by the galaxy oscillation damping time but 
by the tide amplitude (e). For a logarithmic MW potential in a circular orbit the fractional 
resonant-growth frequency domain is simply uj c 2 /ujq 2 (the square of the ratio of the circular 
orbit frequency and the internal galaxy oscillation frequency). Equation 26 implies that any 
dS with oscillation frequency less than about 1.4w c will be tidally disrupted. 

As we noted above, elliptical orbits can be even more likely to disrupt dS stellar systems 
because of their rich epicyclic harmonic structure. Eccentric orbits lead to disruption of 
"stiffer" dS with internal frequencies as large as 1.7, 2.4, and 3.1 (or larger) times the circular 
orbital frequency - corresponding to 1st, 2nd, and 3rd order terms in orbital eccentricity. 
Note that parametric resonance increases the mean dS central distance of its member stars 
and their velocities (and velocity dispersions). The exponential growth time depends on 
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how close the dS is to a resonance condition but as eq. (25) shows the growth times of, for 
example, the velocity dispersion can be comparable to the orbit period. In general we expect 
larger growth rates as the ratio of circular to internal dS frequency increases toward unity. 

Despite the undetermined orbital characteristics of the dS, it is interesting to compare 
our best estimate of e for each of the MW dwarf spheroidals with their kinematic properties. 
It is our contention that the dS kinematics are not dominated by dark matter halos, but 
by their MW tidal interaction. Thus we estimate their mass from their luminosity and not 
their velocity dispersion, since the dS are not in virial equilibrium. Core (half-light) radii 
(r c ) and galactocentric distances (d) are used from the compilation by Mateo (1998). A 
pulsation frequency ojq is computed from the VanderVoort approximation c^o = ^V^Po- We 
compute the central density p from the observed total luminosity (assumed equal to the 
total dS mass in solar units) and the expression p = 3M (°^ . This is exact when the density 
distribution equals a Plummer model and is a reasonable empirical approximation for the 
MW dS. We also take u c = v c /d (with v c = 220 km.s -1 ) as the circular frequency around the 
MW. Adequate for our purposes, we note that the parametric growth rate increases with e 
and the ratio A = uo c /uoq. If the dynamics of the dS are dominated by tidal interactions we 
expect observed dS virial "M/L" ratios to increase with A. Figure 15 plots A against M/L. 
The general trend of increasing "M/L" with A is in good agreement with the parametric 
model, especially given the uncertainty in the dS orbital parameters. 

Our "snapshot" knowledge of the Local Group Dwarf Spheroidals means that we cannot 
measure their orbital ellipticity (although see Kuhn 1993 for a discussion of this point) 
nevertheless, knowledge of the dS radial velocities and galactocentric distances is sufficient 
to suggest that at least the nearest dS must be tidally inflated by parametric resonance. 

Another interesting and perhaps observationally verifiable consequence of parametric 
excitation is that the long axis of the dS bar should, in general, be inclined with respect to 
the separation vector between the MW center and the dS. Orbits which are nearly circular 
should yield a dS with the leading edge closer to the MW than the trailing edge and inclined 
at an angle of 45 degrees. Elliptical orbits produce a rotating bar whose phase with respect 
to the radius vector direction depends on how close the dS is to resonance and the proximity 
of the dS in its orbit to perigalacticon. 

5. Summary and Conclusions 

We have demonstrated an analytic model of time dependent dS-MW tidal interactions 
which extends the Mathieu Equation for parametric oscillations to two dimensions. Several 
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important general conclusions follow from this model: 

1. Secular instability can cause dS star orbits and velocities to grow, depending on the 
ellipticity and circular frequency of the MW orbit and the resonant frequency of the 
dS. 

2. Exponential growth can occur over a broad range of frequencies, not just near the 
characteristic dS harmonic frequency. Growth times can be comparable to the MW 
orbital period. 

3. Direct integration of more realistic, non-harmonic dS potentials, yield the same quali- 
tative behavior predicted by the Mathieu equation. 

4. Self-gravitating N-body calculations also confirm the existence of broad resonance con- 
ditions that lead to stellar ejection from the dS system, elliptical bar formation, and 
inflated (non-virial) velocity dispersions. 

5. Model expectations for the parametric growth rate in MW dwarf galaxies have been 
contrasted with the observed M/L measurements for the 8 Milky Way dS. The expected 
trend of increasing M/L with u c /uo is confirmed. 

Considerable attention has been given to finding dark matter density distributions which 
will both stabilize the luminous components of the dS so they survive MW tidal encounters, 
but not so much that they don't allow stars to be ejected to populate "tidal tails" . These 
tails are now clearly observed by several methods (Kuhn et al. 1996, Smith et al. 1997, 
Martinez-Delgado et al 2001) and are an observational hurdle that any dS model must pass. 
While it appears that there may be dark matter models which cannot be ruled out (cf 
Mayer et al. 2002) here we have developed a consistent model of MW-dS interactions which 
does not require any invisible dS mass component and which appears to account for the 
morphology and dynamics of the dS. The prediction of extratidal dS stars in the simple KM 
resonance model were later observed. We can hope that the considerably refined predictions 
of the parametric oscillator model (for example the inclined bar) may also be empirically 
established. 
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Fig. 1. — Tidal interaction geometry A dS star is represented by 8r in the dS-frame 
centered on C and by f + Sr in the MW-frame of reference centered on O. 
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Fig. 2. — 6 = oo, harmonic case . Decoupling the tide amplitude and circular frequency 
allows us to see where the parametric excitation takes place. For this Runge-Kutta integra- 
tion over 50 time units, we used ojq = 4, b = oo, e = 0.1 with initial conditions x(0) = 1 
and v — 0. The resulting RMS stellar orbit radius peaks for dS circular frequencies near uj 
with a half-width consistent with euj /2 = 0.2 as predicted by the analytical solution. Note 
also that the one-dimensional parametric oscillator excitations at uo/n with n > 2 do not 
appear, as predicted by our analytic model. 
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Fig. 3.- 



0, unbound case. Symbolic computer aided analytic solutions for this 



unbound case suggest that the dominant term should be of the form exp \/eu>o 2 — uP t, which 
gives a cut-off frequency u cnt = ^/eooq = 1.26 as observed here. Stellar orbits are stabilized 
(bounded) when uj > y/eujQ = y/k where k is the tidal constant defined in eq. (11a). 
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Fig. 4. — r/b = 1 (top panel) and r/b = 1.4, (lower panel) intermediate cases. 

Conditions are identical to figure 2 except that b has now be chosen so that x(0)/b = 1 and 
x(0)/b = 1.4. The apparent excitation frequency decreases as you move further away 
from the core. 
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Fig. 5. — Elliptical orbit excitation. Runge-Kutta integrations for a dS with c^o = 6.7, 
x(0)/b = 1.4 in elliptical orbits with fixed perigalacticon distance. Here e = 1 — u/2 so that 
uj = 2.0 is a circular orbit. The tide amplitude and circular frequency are constrained as 
they are in a logarithmic potential so that k = uj 2 . Orbital and dS parameters correspond 
approximately to the N-body simulation parameters. The broad frequency response of the 
dS is a consequence of elliptical orbits. 
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Fig. 6. — Influence of the softening parameter. The potential softening parameter also 
affects the relaxation time and numerical dissipation. Here we plot the galaxy oscillation 
spectrum for identical initial N-body conditions but with softening parameter values or 
a = 0.025 (upper panel) and a = 0.1 (lower panel). Curves are generated from the timeseries 
defined in equation (34). 
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Fig. 7. — Influence of integration time-step on relaxation. To produce long simula- 
tions where numerical dissipation is minimized, the integration time-step must be decreased. 
These curves show the evolution of isolated sytem rms radius for 3 runs with integration 
time-steps of 1/32, 1/128 and 1/1024 time units. The radius of the system grows apprecia- 
bly for 1/32 time-steps but is essentially stable for the smaller integration times. We adopt 
1/128 in our integrations but note that the fine structure in these calculations (consider the 
difference between the two shorter integration rms curves) may not be reliable. 
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Fig. 8. — Snapshots at uj p /2 (upper left), 2uj p (upper right) and uj p (lower panel). 

Snapshots of the projected dS star distributions have been taken at t = 25 numeric units 
for three runs (uj = 0.8 (upper left), uj = 4 (upper right) and uj = 1.7 (lower panel)) 
corresponding to uj uj p /2, uj ks 2uj p and uj ks uj p . The instantaneous direction of the tide 
is indicated by the solid line. Snapshots at uj p /2 and 2uj p are very similar: the dS does 
not seem to be affected by the tide. The apparent radius is smaller for case uj p /2 because 
almost a third of the stars have been ejected to infinity, although even in this case, the 
dS shape is not affected. On the contrary, the snapshot at uj p shows a dS that has been 
completely disrupted by the tidal interaction. Its elliptical shape is turning with the tide, 
but keeping an approximately constant angle with the line-of-sight - consistent with the 
parametric predictions for circular orbits. 
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Fig. 9. — Number of particles lost versus circular orbit frequency. These simulations 
were performed using a simple turning tide with a constant amplitude k = 2 corresponding 
to e 0.5. Particles are defined to be lost when they reach a distance greater than 1 kpc 
from the center of mass position. The multiple curves plot the number of particles lost 
at successive times. The curves are separated by 6.25 time units. The effective resonant 
frequency of the dS is uj v = 1.7 initially. As stars are lost and the central density decreases, 
uj p decreases and we also observe the resonant peak shift to lower frequencies. We also 
observe two regimes in the simulation results - if to < 1.5, lost stars are ejected to reach 
r = oo, but when oj > 1.6, stars remain in a bounded region of space near the dS. 
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Fig. 10. — Parametric excitation on elliptical orbits. This figure plots the variation 
in the number of particles lost versus time and orbit frequency. Multiple curves in the figure 
show the dS system at successively later times At = 25. With a fixed perigalacticon distance 
at 50 kpc we have u c ~ 2 (1 — e), so that oo c = 2 corresponds to the circular case e = 0. The 
maximum loss is 1024 stars but all simulations began with 734 stars when injected into their 
MW orbit. Elliptical orbit tidal interactions can expand and entirely disrupt a galaxy due 
to parametric excitation where a stronger static tide has only a small effect. These results 
compare favorably with the non-gravitating Runge-Kutta elliptical orbit calculations above 
(c/ figure 5). 
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Fig. 11. — Velocity dispersion. This figure plots the derived line-of-sight velocity disper- 
sion near a weakly resonant system with orbital eccentricity 0.5. The fluctuating velocity 
dispersion can be many times larger than its initial equilibrium value before the dS appears 
to be disrupted. Virial mass estimates near the perigalacticon dS passages could yield M/L 
values exaggerated by a factor of 100. 
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Fig. 12. — Snapshot near a maximum velocity dispersion The dS stellar distribution 
as projected against the sky is shown here. The plane of the orbit is indicated by the solid 
line. This view corresponds to the dS when the velocity dispersion was approximately 10 
times larger than its initial value and the core of the dS is apparent even with a small number 
of stars in our N-body calculations. 




Fig. 13. — Snapshots of major axis rotation. Snapshots of the dS stars projected 
into the orbital plane are shown here. Time advances from left to right and top to bottom 
with each panel 0.25 time units apart. The horizontal line represents the direction of the 
line-of-sight and the solid line shows the direction of the major axis as computed from the 
moment of inertia tensor of the stellar distribution. Rotation occurs from snapshots 3 to 9. 
Note also that stellar ejection occurs primarily when the major axis begins to turn, ejecting 
groups of stars into the tidal tails. 
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Fig. 14. — Variation of the angle at first order in elliptical excitation. Rotation of 
the dS bar appears as a discontinuity in the derived bar angle. We have not represented the 
first 50 time units where the major axis is not well defined and the angle varies chaotically. 
Between successive perigalactica we see one discontinuity corresponding to half a turn of the 
dS. This is characteristic of excitation at first order in e. 
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Fig. 15. — Mass/Light ratio against estimated growth rate. The Local Group dwarf 
spheroidal galaxy virial "Mass/Light" ratio in solar units is plotted against a derived mode 
growth rate (A = uo c /uoq). Parametric oscillations should inflate the velocity dispersion (or 
M/L) with increasing A as is seen here. 



